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Abstract 

> 

^ I We derive a system of covariant single-time equations for a two-body bound 

^ ! state in a model of scalar fields 0i and 02 interacting via exchange of another scalar 

lO ' field X- The derivation of the system of equations follows from the Haag expansion. 

masses, mi and m2, of the scalar fields, 0i and 02- We present an approximate 

Oh: 



The equations are linear integral equations that are explicitly symmetric in the 



analytic formula for the mass eigenvalue of the ground state and give numerical 
results for the amplitudes for a choice of constituent and exchanged particle masses. 
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1. INTRODUCTION 

The problem of relativistic bound states has a long history. Nonetheless, the 
treatment of this problem is still not completely satisfactory. The purpose of this 
paper is to continue the development of an alternative to the most popular formu- 
lation, the Bethe-Salpeter method[|l[. The Bethe-Salpeter method uses amplitudes 
in which both constituents are off-shell. Because of this, the amplitudes depend on 
an unphysical relative-time coordinate and obey equations that are difficult to solve 
and have spurious unphysical solutions, including some of negative norm. Several 
authors have proposed covariant, single-time equations with only one constituent 
off-shell. The equations most similar to ours are the "spectator equations" of F. 
Gross|0]. Our equations differ from the spectator equations in the way we ensure 
symmetry between the off-shell and on-shell particles, in the inclusion of renormal- 
ization graphs and counter terms and in the boundary conditions of the Green's 
functions: our equations use Green's functions with retarded boundary conditions, 
rather than Feynman boundary conditions. Our derivation of the equations dif- 
fers entirely from Gross' derivation of the spectator equation: we use the Haag 
expansion[^ and the operator field equations^, ^, |^, H, rather than summing 
classes of Feynman graphs. To be concrete, we consider a two-body bound state 
in a model of scalar fields 0i and 02 interacting via exchange of another scalar x- 
A. Raychaudhuri used this method to study bound states in the equal-mass case of 
this model ; however his equations are not symmetric in the on-shell and off-shell 
masses0. This asymmetry was not evident in the equal-mass case. He also stud- 
ied the nonrelativistic reduction of the equations for unequal-mass bound states of 
spin-1/2 particles|T^. Related work was done by M. Bander, et al|ll]. 

In this paper, we extend Raychaudhuri 's analysis to unequal-mass constituents 
and treat the on-shell and off-shell particles in a completely symmetric way. We 
present numerical results for the ground state eigenvalue and amplitudes for a range 
of mi/ 1712 and yu/m2, where /i is the mass of the x field. 

We hope this method will provide an alternative to the Bethe-Salpeter method 
with the following advantages: (1) because only one particle at a time is off-shell, the 
amplitudes depend on only one invariant, (2) all normalizable solutions are physical 
and have positive norm, (3) the limit for one mass very large is the relativistic 
equation for the other particle bound in an external field and (4) the nonrelativistic 
limit has the correct reduced mass. We have two longer-range goals: (1) to extend 
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the Haag expansion to account for cases in which the interaction is strong enough to 
make a two-particle description inadequate and (2) to modify the Haag expansion to 
treat confined degrees of freedom, for which the usual asymptotic fields don't exist. 
This latter goal will require a significant generalization of the method. 

2. DERIVATION OF THE EQUATIONS 

In this section we obtain coupled integral equations for two bound-state am- 
plitudes, one, /i, with particle one off-shell and particle two on-shell, the other, /2, 
with these roles interchanged. Our Lagrangian is 

^ = 12 lid,<P,d'^<P,) - ml<Pl) + \{d,xd''x - /^V) + f [0? + 4, X]+, (1) 

where the last term is an anticommutator. We work in momentum space using 
= (27r)~^/^ / d^p^{p)exp{—ip ■ x) and the analogous formula with in fields. We 
promptly drop the tilde on 0, abbreviate d^p by dp and, for the in field, abbreviate 
■ 4>^"'{p)^{p^ ~ mf) : by : :. The equations of motion in momentum space are 

{m] - p^)(t)i{p) = ^^^^ j dpidp25{p-pi-p2)[(l)i{pi),x{P2)]+ (2) 

- Biml)(t>i{p) (3) 

g 2 

(/^^-P^)x(p) = J dpidp2S{p - Pi - P2)(l)i{pi)(f)i{p2) (4) 

+{Dp' - Ef,')x{p), (5) 

where we have introduced counter terms for the mass and field strength renormal- 
izations of the fields. 

In the A'"-quantum approximation, we expand the Lagrangian fields in terms 
of the complete, irreducible set of in-fields (or out- fields), including those for stable 
bound states (this is the Haag expansion), and truncate the expansion to find an 
approximate set of equations among a finite number of amplitudes. Here we keep all 
terms that contribute to equations for the two-body bound-state amplitudes fi and 
/2 mentioned above in one-loop approximation. All the terms have explicit order 
at the perturbative vertices. The relevant terms in the Haag expansion are 
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Mp) = : ^tiP) ■■ + I dqdbdip + q - b)Mq, h) : (t>^{-q)B'-{h) : 

+ j dqdldh5{h -q-l- h)f^2B{q, h, I) : x''^(g)C(0^'"(^') 

xip) = ■.x'''{p)-.+ j dhdh5{p-h-hHi{hM)--^nhW^^^ 

+ j dhdkdbSip -h-h- b^uBih, k, b) : <f)f{h)(f>T{l2)B'''{b) :, 

and the analogous terms for 02, with 1 and 2 interchanged. Our general notation is 
that the subscript on an amplitude lists its associated product of in fields. What we 
call fi should be according to this general notation; however, for convenience, 
we call it /i. Because each in field has a mass shell 5-function, all the momentum 
integrals in the Haag expansion are on two-sheeted mass hyperboloids. In fi we 
keep b on the positive-energy mass shell and reverse the sign of the momentum q 
in /i so that q on the positive-energy hyperboloid gives the dominant amplitude in 
the nonrelativistic limit. We call fi{q,b) with q > 0, i.e. with q on the positive 
mass hyperboloid, f^^^ and with q < 0, fi~\ (See Fig. 1) Both our equations and 
their graphical representation include both of these pieces of the amplitudes; to save 
space we don't exhibit both pieces in the graphs. 

As usual, : : denotes normal ordering. In the one-loop approximation, contrac- 
tions always involve the vacuum matrix element of the anticommutator, 

([C(Pl)>f(^^2)] + )o = S,j5(p,+p2)5mXPl), 

([x^"(pi),X^"(p2)]+>o = Sip,+p2)5M, 

where Sm{p) — S{p'^ — rn^) for short. Choosing to expand the Lagrangian fields in 
terms of the in-fields requires using retarded boundary conditions for the propaga- 
tors. 

To obtain the equation for /i, we insert the Haag expansions for 0i and x 
in the equation of motion for 0i, renormal order and equate the coefficients of the 
term with : (f)^^B^"^ :. The resulting equation involves the amplitudes /i, fx2B, I22, 
and 7i2B. We calculate the last three amplitudes in terms of /i using the equations 
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of motion and the Born approximation for emission of both on-shell and off-shell x 
quanta. We will give details of this in a later, more detailed paper. 



/x2B(g,^,&) 

7ll(Pl,P2) 
ll2B{hj2,b) = 

The integral equation for fi is 
[ml-{b-p)']f,{p,b) 
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fi{-l,b) 



(27r)3/2 ml-{q + l + b)'^' 
9 



722(Pl,P2 
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2(27r)3/2 f^2_^p^^p^y 

fii-l2,b) + U-h,b) 



(6) 
(7) 



(27r)3/2 ^2 _ (/^ + + 



(9) 



9' 



167r3 



dq 



+ 



U1) 



IJ? — [b — p — qY rrii — {b — p — q) 



167r3 



dq 



(p - q) 



;fiiq,b) + 



(b-p-q) 



fi{q,b) 
f2{q,b) 



+ [A,ib-pf-B,ml]Mp,b), 



(10) 



where the first two terms on the right hand side are self-energy graphs that are com- 
pletely canceled by the renormalization counter terms. Any method of regularization 
of the self-energy graphs will suffice. The third and fourth terms on the right give 
binding by exchange of the x field. The bound-state momentum b is always on its 
mass shell 6^ = M^. This bound-state equation is shown in Fig. 2. We get another 
coupled equation for /2 by interchanging 1 and 2. The resulting pair of equations is 
clearly sjTiimetric under 1, 2 interchange. We suppress the ie's associated with the 
retarded boundary conditions. 

3. Approximate Mass Eigenvalue Formula 

We considered parametrizing the mass eigenvalue formula using the arccosT], 
where 1] = M/{mi + 1712), because this expression appears in the hyperboloidal 
harmonic analysis that Raychaudhuri]^] used. We found interesting empirical regu- 
larities using this parametrization. The result is 

A — a 



M = (mi -|- 1712) cos ■ 
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where A = (?^/(327rmim2), a = Q.9\J^/mred and b = 0.8 — 1.1 ln(m</m>). The 
reduced mass is the usual expression; m> is the larger of mi and 1712. The range of 
validity of this empirical formula is < /i < m<, 0.01 < m^/rriy < 1, 0.5 < r] < 1 
for m^/niy = 1 and 0.9 < rj < 1 for m^/m^ = 0.1. 

4. Numerical Results 

Equation (9) and the one with mi and m2 interchanged are eigenvalue equa- 
tions for the coupling constant g; that means for given values of the masses M, mi, m2 
and /i we can find a coupling constant g and wave functions /i^^ that satisfy the 
equations. We solve these homogeneous linear integral equations by approximating 
the integral on the right hand side with a finite sum. We choose Gauss integration 
with appropriate points and weights. The resulting matrix equation is solved by 
standard means. 

For the equation in momentum space it is sufficient to take 18 mesh points to 
obtain g^ to an accuracy of 4%. The main difficulties we encounter in Eq. (9) are 
the logarithmic singularities. We smooth these singularities by keeping a finite e at 
the logarithmic singularity. We checked that the result does not change by varying 
the mesh points and e. 

In Fig. 3 the value of A = g"^/ (327rmim2) is plotted as a function of rj = 
M/(mi +m2) for m^/m^ = 0.1 and /i = 0. A calculation using the Bethe -Salpeter 
equation |]l2l consistently gives smaller binding. In the scalar model we cannot 
decide which solution is correct because there is no experimental data. Figure 4 
shows the wave functions fi^\ fi^\ /i^'' and f2~\ respectively, for the mass ratios 
given above and for rj = 0.95. As expected /i^"* is the dominant contribution. 

We will present more extensive numerical results in a later paper. 

5. Summary and outlook for future work 

The Haag expansion leads directly to coupled linear integral equations for four 
amplitudes related to a scalar bound state. Two amplitudes are those that reduce to 
the nonrelativistic wavefunction, one, /i, with the particle of mass mi off-shell and 
the other, /2, with the particle of mass m2 off-shell. The other two amplitudes have 
the on-shell particle crossed, so that its momentum lies in the same light cone as the 
bound-state momentum. These four amplitudes obey a set of four coupled linear 
integral equations. We solved these numerically using momentum-space variables. 



6 



We plan to apply this method to bound states of two spin- 1/2 particles, such 
as the hydrogen atom and positronium, where our calculations can be compared 
with experimental results. We hope this method can replace the Bethe-Salpeter 
method in theories without confinement. 

In order to use this method in confining theories, such as QCD, the asymptotic 
fields that are a prominent part of the Haag expansion must be replaced with fields 
that correspond to confined degrees of freedom. The treatment of confined degrees 
of freedom in theories such as QCD remains a goal for the future. 
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Captions 

Fig. 1: The two pieces of /i. The short hne through a leg indicates the leg is 



Fig. 2: Graphs for the bound state equation for fi if the left-hand leg is 0i. 
Note that the first two terms are self-energy graphs that are canceled by the counter 
terms, the third term is a i-channel graph that couples /i to itself and the last term 
is a li-channel graph that couples /i to f2- 

Fig. 3: Plot of A = g'^ / {327rmim2) as a function of 77 = M/{mi + 1712) for 
H = 0, m<:/my = 0.1. 

Fig. 4: Wave functions fi^\ f[ \ and /2 "* in momentum space in ar- 
bitrary units as a function of A for 77 = 0.95, = and m</m> = 0.1. Here 



off-shell. 
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Fig. 3 
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Fig. 4 



